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It is shown that hyperspherical harmonics can be represented in a form typical of traditional 
microscopic cluster models. This allows those hyperharmonics that are responsible for long-range 
behaviour of valence nucleons in loosely-bound nuclei to be selected. The hyperspherical cluster 
model based on such hyperharmonics is tested for ^He using a simplified description of the *He core 
and two versions of the Volkov potential which give ^He as either bound or unbound. The study has 
confirmed that it is possible to get a converged solution for binding energy, r.m.s. radius, overlap 
integral and asymptotic normalization coefficient with small number of basis functions for bound 
^He. The values obtained of these quantities are very close to those calculated in a microscopic 
cluster model. 
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. Weak binding in nuclei with a large excess of neutrons or protons leads to strong clustering and to a slow decrease 
^Nj ■ of the wave function of valence nucleons at large distances. To describe the processes in which such nuclei are involve 
04 I and to deduce structure characteristics through the study of these processes, long-range behaviour of valence nucleons 
should be properly taken into account. Current reaction theories often use simplified phenomenological nucleon-core 
" ' ' models to include long-range behaviour relevant to reaction mechanisms. However, a deeper understanding of nuclear 
^ I ■ structure and its manifestation in nuclear reactions requires nuclei to be considered as many-body objects. 
\^ I It is not easy nor always possible to obtain correct long-range behaviour from first principles. It can be obtained 
. naturally only for three-body and four-body systems using Faddeev and Faddeev-Yakubovsky equations. For three- 
T-H nucleon systems, correct long-range behaviour has also been achieved using pair correlated hyperspherical harmonics 
. which has allowed asymptotic normalization constants for the n — d and p — d systems to be calculated. The 

hyperspherical formalism has also been applied to calculations of the overlap integral between triton and ''He 0j. 
Despite the very large basis space used in these calculations, the proper asymptotic behaviour of the overlap (tig)p|'*He), 
given by the Whittaker function with relevant proton separation energy, has not been achieved at r > 6 fm. To restore 
this behaviour, the authors solved the inhomogeneous equation for overlap integrals fro Ref. For heavier nuclei, 
only a few ab-initio calculations of overlap integrals are available. Large-scale no-core shell model calculations of the 
overlap (^Be®p|^B) [2| have not provided the required decrease at r > 5 fm so that the authors had to match this 
(~| , overlap to the Whittaker function, or approximate it by solutions of a Woods-Saxon potential in order to use this 
overlap in reaction calculations. Three other ab-initio overlap calculations are available within the variational Monte- 
Carlo method, for {'^Re(g)dfU) C^He^tl^Li) and (^He^^Hel^Be) 0. The radial behaviour of the Ion g-range part 
of these overlaps is imposed there prior to calculations so that it is not determined by the differences in calculated 
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5h , energies of the nucleus and the clusters that compose it. 

. . 1 Many-body approaches that use Slater determinants made of single-particle wave functions have also to adjust the 
potential wells to fit the experimental energies for valence nucleons if such wave functions are to be used in reaction 
calculations. However, this procedure does not necessarily provide the separation energy that is the difference between 
the total energies. Also, it may be inconsistent with nucleon-nucleon (NN) interactions employed in these approaches. 

Proper long-range behaviour can be guaranteed if the total wave function is represented by an antisymmetried 
product of the core 0c and valence (j^yai wave functions 

* = A{(f>c ® (pval) (1) 

and then the microscopic R-matrix approach is used determine the function (f>vai [ll- present a microscopic 
cluster model (MCM) of this type is used to predict cross sections of astrophysically relevant low-energy reactions [8] . 
However, this model is only able to deal with oscillator shell model core wave functions (t>c and can not be used with 
realistic NN interactions that reproduce NN data. Although in this model the separation energies are always equal 
to the differences of the total binding energies, they often differ significantly from experimental values so that one 
parameter of the effective NN potential should be tuned. 

The problem of the long-range behaviour can be resolved within the hyperspherical interpolation approach proposed 
in [P]. This approach has been designed to describe many-nucleon systems near thresholds when identified channels 
having a binary cluster structure are believed to play an important role, especially in the nuclear surface and beyond. 
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In this approach the wave function of a many-nucleon system ^' is represented by two terms ^' = + \E'2, one of 
which exphcitly contains the binary channel wave function. Both and '^2 are expanded in hyperspherical basis and 
a system of equations that couples them derived by substituting the total wave function into the general expression 
for the variational principle. The interpolation approach does not involve matching at some surface. The coupling 
between the two components of the wavcfunction is governed by the NN interaction through the Schrodinger equation. 
Earlier applications of the interpolation method included the calculation of ^He+n elastic scattering [ifl]- Later, this 
approach was used to study the long-range behaviour in *B and ^Li nuclei within a three-body model S imilar 
ideas have been successfully used to describe the '^He-l-p scattering within the hyperspherical formalism The 
success of these approaches suggests that the hyperspherical interpolation method can be applied to heavier systems 
as well. However, to make it applicable to other systems, a hyperspherical expansion that contains contributions only 
from the long-range part of the wave function ^E" should be introduced. 

In this paper, a subset of hyperspherical basis is constructed for A identical fermions that is expected to be 
responsible for long-range behaviour in weakly-bound nuclei and a technique to calculate matrix elements in such 
a basis is developed. This basis is used to expand the total wave function 4" into hyperradial and hyperangular 
parts and then the hyperradial part is found by solving the standard set of coupled differential equations of the 
hyperspherical functions method (HSFM) . The feasibility of the proposed method is studied in the simplest case of 
^He. Only those hyperspherical basis states have been chosen that correspond to the ^He core described in the lowest 
order approximation of the HSFM. The convergence of such an expansion is investigated by calculating the binding 
energy and r.m.s. radius of '^He with an increasing number of basis states using the Volkov effective NN potential 
VI [161 with two values of the Majorana parameter, the standard m — 0.6 that corresponds to unbound ^He and the 
non-standard m = 0.3 that binds this nucleus. For the latter value, the convergence of the hyperspherical expansion is 
also studied for the overlap integral (^He(8)n|^He), its r.m.s. radius, spectroscopic factor and asymptotic normalization 
coefficient (ANC). The hyperspherical expansion of the present paper should produce results that are very close to the 
MCM with a closed Os-shell "'He core described in the oscillator shell model. The comparison between both models 
is presented. 

In Sec. II the hyperspherical cluster basis is constructed and in Sec. HI the expansion onto this basis is discussed. 
The link between the hyperspherical cluster model and MCM is shown in Sec. IV, while the representation of 
hyperspherical cluster harmonics by oscillator cluster wave functions is given in Sec. V. A method to calculate matrix 
elements in this basis is presented in Sec. VI and the application to ^He is discussed in Sec. VII. The results obtained 
are summarised in Sec. VIII. Finaly, exact formulae for the norm and the two-body NN potential matrix elements as 
well as for the overlap integral are derived in the Appendix. 



II. HYPERSPHERICAL CLUSTER HARMONICS 



A nucleus A can be described in hyperspherical coordinates, that include the hyperradius p, 

p^-iie^ (2) 

where are the normalised Jacobi coordinates, 

and a set of 3A — 4 hyperangles p. The wave function 5* written in these coordinates can be expanded into hyper- 
spherical harmonics (HHs) that depend on the hyperangles. 

Let us consider a nucleus A that is strongly clusterized as a core A—1 plus a weekly-bound valence nucleon N. It is 
well-known [l^ that the spatial HHs for A nucleons can be written as a product of a spatial HHs for A—1 nucleons, 
the spherical function ilm(^A-i) ^ function ipnid) which is an eigenfunction of the angular part of Laplacian in 
variable 9 defined as 6* = arctan^^_i/pc- Here pc is the hyperradius for the A—1 core, = p1+ The angle 9 

is associated with the valence nucleon and it shows how far the valence nucleon is for a fixed size of the core A—1. 
To describe the long-range behaviour of the last nucleon for a fixed structure of the core it is necessary to include as 
many eigenfunctions ipniO) as possible. 

This paper accounts for cluster structure (yl — 1) + A'^ of weakly bound nuclei by reorganizing known recurrence 
representation of HHs in a form that is similar to the microscopic cluster model. This is achieved by antisymmetrizing 
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the product of a (known) completely antisymmetric HH y^^^c^^^c^^^c ^p^^j fQj- ^j^g core A^l (spin-isospin part included) 
and a relative angular function ^K^nimaT{OiS,A-i) 

'>^KaianlLST\r) ^ ■'^Kai^nlLST / . 

x{L,MLjm\LML){ScMs^^a\SMs){T,MT^^T\TMT)A (y^'X^"'^^'^^ {p^) VK^nlm.r{0,iA-l)) ■ (4) 

Here is the hypermoment, 7c — {PcLcScTc}, Lc, Sc, and Tc are the total orbital momentum, spin and isospin of 
the core A — 1 while the /3c represents different degenerate harmonics. Also, 

is the antisymmetrisation operator which permutes the A-th nucleon with nucleons of the core, MK^-y^ni lst is a 
normalization factor, denotes the hyperangles for the A ~ 1-body system and 

<PK^nlmar{e,kA-l) = NnKJ (siu^)' (cOS0)^= P^+l/2.^= + (3^-8)/2 (^^^ 20)y,™ (^^_ 1 )X.r ( A) (6) 

where cos0 = Pc/p, sin 6 = S,a-i/p, cos 20 = 1 — 2^\_^l p^, X<tt{A) is the spin-isospin function of the a-th nucleon 
with the spin and isospin projection a and t respectively, and 

2 ^ 2n! (2n + + ? + (3^ - 5)/2)r(n + K^ + l + jiA- 5)/2) 

r(7i + Z + 3/2)r(n + /fc + (3A-6)/2) ' ^' 

Since Y^^'^"^^'^"^^^" (p^) and ^Kanima-T{0,^j^_i) are the eigenfunction of the operator of the kinetic energy in sub- 
spaces associated with variables p^ and {^a-it^} ^^'^ because the operator of the kinetic energy is symmetric and 
the function ^ is antisymmetric with respect to any nucleon permutations, the function yK^^^niLST^P) = yK-yip) 
is also an eigenfunction of the 3 A — 4-dimensional angular part Ap of the Laplacian, 

ApyK^{p) = -K{K + 3A~5)yKy{p), (8) 

and therefore it is the HH with the hypermoment K ^ Kc + 2n + l. Here 7 = {Kcjcnl LSTMlMsMt}- 

Since the HHs yx-y (p) are written in the form typical for the cluster wave functions ([ij , they are referred to below 
as hyperspherical cluster harmonics (HCH). The HCHs yK-y{p) with different values of K are orthogonal to each 
other, however, for the same K and different 7 they may be not orthogonal 

{yK'Y{p)\yK'y{p)) ^ SkK'^KjY, (9) 

which is common for different channel functions in multichannel cluster models with antisymmctrization. 



III. EXPANSION OF NUCLEAR WAVE FUNCTIONS ONTO HYPERSPHERICAL CLUSTER BASIS 

The wave function of a nucleus with one weakly-bound nucleon can be expanded onto the hyperspherical cluster 
basis, 

^-p-^'^-'^^'Y.^K,{p)yK,ip), (10) 

keeping only a restricted number of quantum numbers in the HHs of the core and as much relative functions ([6]) (or 
in other words as many n's) as needed to describe properly the long-range radial behaviour of the valence nucleon. 
However, when more then one set of quantum numbers for the HH of the core is present in the expansion (I10|l . the 
hyperspherical cluster basis may not be orthogonal for the same K. The simplest way to deal with this problem is to 
introduce the orthogonalised HHs, 
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(11) 
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which can be achieved, for example, by performing the singular value decomposition of the matrix XK-y-y' ■ Then the 
decomposition 

= p-^'^-'"'Y.^K,{p)yK,{p) (12) 

Ky 

onto orthogonal basis yxjip) leads to a standard set of differential hyperradial equations of the HSFM, 
(P £k('Ck + 1) 2m -- / nn\ ~ / N 2m ^ ~ 

-p^ 2 -rri^ + ^K^>K^[p)) j XK^[p) = ^ 2^ Vk^,k'^'{p)xk'j'(p), (13) 

ap P ti J h j^,^,^^^ 

where Ck = -ftr+(3yl — 6)/2, m is the nucleon mass and the hyperradial potentials Vxy.K'y'ip) are the matrix elements 
of the NN interactions 

VK^,K'y{p) = {yKy{p)\ J2 ^^li^^ - r,)\yK'Ym- (14) 

The functions Xk-i{p) corresponding to the non-orthogonal HCH basis yxjip) with well-defined quantum numbers 
{Kcjci^l} are then obtained using the transformation 
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IV. LINK BETWEEN THE HYPERSPHERICAL CLUSTER MODEL AND MICROSCOPIC CLUSTER 

MODEL 

If expansion (llOp retains only a limited number of {Kc'jc} while at each fixed {Kcjc} the summation over all 
quantum numbers n is performed, such an expansion corresponds to a specific cluster model of the type ([T]). The 
easiest way to demonstrate this, is to retain only one state {Kcjc} for the core A — 1 and neglect all spin couplings. 
Then the wave function (IIOII can be rewritten as 

(16) 



* --^ i^^K^^j^nlLST ^J3A-l)/2 ^fcc7c(Pc) 'fiK^nlmariO , 
\ n ' 



If after the antisymmetrization operator in (|16|) a closure relation is inserted 



KUSPc,Pc)){KUSPc,Pe)\^l, (17) 



where ^'^''^ (pcPc) is the i-th solution of the Schrodinger equation for A — 1 nucleons in the hyperspherical basis 
that keeps only one selected HH, Yk^^^{p^), then the wave function PH)) can be rewritten as 

* = J2^{'^kL(Pc,Pc) 4[,Ma-i)) , (18) 

i 

with the relative function 4>\^ ^ (€a-i) determined by the following expression: 

Thus the HCH expansion of the wave function ^! is equivalent to a microscopic multichannel cluster model that 
includes all the excited states of the core A—1 that correspond to the hyperradial excitations of the lowest state with 

{ifc7c}. 

In the general case, when more sets of the quantum numbers {Kc'^c} are retained in the HCH expansion, it is 
possible to show that the HCH expansion corresponds to the multichannel cluster model that includes all the excited 
states of the core with total spin J^, that can be constructed with the values of and Sc retained in the core 
description, and their hyperradial excitations. Such a derivation is quite cumbersome and is not given here. 
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V. LINK BETWEEN HYPERSPHERICAL CLUSTER HARMONICS AND OSCILLATOR CLUSTER 

WAVE FUNCTIONS 

The HCH yK-y{p) can be represented by the cluster wave functions that contain oscillator shell model wave functions. 
Such a representation is used below to develop a technique to calculated the matrix elements in the HCH basis. The 
link to the cluster shell model can be derived using an expansion of the Jacobi polynomials into the radial oscillator 
wave functions (see Eqs. (|47l) - (|49|) of Appendix A) and representing the core HH YK^-y^iPc) as 

YK^yApc) - '^okM^i,-,^a-2)/RokAPc), (20) 

where ^'oifc7c(^ii ■■■1^71-2) translation-invariant oscillator shell model function without the hyperradial excita- 

tions, 



is the hyperradial oscillator wave function and b is an arbitrary oscillator radius. For the lowest possible value 
of Kc the ^'oA'e7c (^iJ ^^-2) ^^e usual Ohuj translation-invariant shell model wave function. For higher Kc, 
^0Kc7c (^iJ ^A-2) can be constructed, for example, using the technique of Ref. [14|. With Eqs. pT|) . ((20l) and 
dSJ-dlSl) one gets 

n 

yK^ip) = AA^; E E ^nKJ ZaAp) (22) 

U = K=0 

where a = {KcIchILSTMlMsMt}, 
and ^'i\ii,..4A- 

is a non-antisymmetrized oscillator cluster wave function, 
^'-'\^i,-4a-i) = E {LaMLjm\LML){ScMs^^a\SMs)iTaMT^^T\TMT) 

X*olf';'/^=*'"^(€l,...,^A-2)i?K;™(^A-l)X.r(A), (24) 

corresponding to the oscillator radius b. The expansion coefficients -B^^ ; are given by the expression 

1/2 



(23) 



^nKJ = (-) J^nKJ ^^^j 

+ / + 1/2 \ / n + ; + 1/2 \ f n + Kc + l + i3A- 7)/2 + z/ 
K + l + l/2)[ n-v )\ n + Kc + l + {'iA-r)/2 



(25) 



VI. MATRIX ELEMENTS IN THE HYPERSPHERICAL CLUSTER BASIS 

The matrix elements in the hyperspherical cluster basis for an arbitrary operator O can be calculated using the 
link between the HCH and the oscillator cluster wave function. This gives 

n n' u u' 

{yK'^'{p)\d\yK^{p)) = {Uk'^,Uk^)-'Y.11 E E ^n'5f^;'^«K./ {Zo^'.'{p)\d\Z^,{p)) 

v=0 i/'=0 K=0 k'=0 

(26) 

To calculate the matrix elements ^^qv (p)|0|^Qiy(p)^ the technique of Ref. fl^ can be used that replaces the 
integration over hyperangles by the Laplace transform of the shell model matrix elements. This results in 

(Z„v'(p)|0|^a.(p)\ =p-^=-'^^-'(''+"'^-'-''-='^+'- / dse^^%"(^=+^'^+'+''+'^-^-3)/2-.-.'oa (^-1/2) (27) 
\ / J -IOC 
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where 

= J di,...d^^_,A ..4a-i)) OA (<i>i^)(^i, ..4a-i)) 

(28) 

are the translation-invariant cluster shell model matrix elements corresponding to the oscillator radius b = s^^. 

A. Normalization of hyperspherical cluster harmonics 

The normalization coefhcients Afx-y of the HCH are obtained using Eqs. ([9]), (I26p - (l28l) with O = 1. Since the norm 
^aa' = {A^a'\A^a) of the osciUator shell model cluster functions does not depend on b the integration over s in Eq. 
(P7)) can be performed immediately. This gives 

^ n n v v' X ryu' k' tjvk tSM 

^ _ 2 \^ \^ <^K'^+2k'+1',K,+2k+1 ^nK'J'^nKJ^aa' ,^9) 

~ AfK-yJ^K-y' ^ ^ ^ ^ p { K,+K'^+l+l'+3A-3 ,\ ' 

j^j j^i j,^oiy'=OK=OK'=o i I 2 \- 1^ + ly 1 

Then AfK-y is found from the condition Xxj-f = 1- 

It is useful here to clarify the meaning of this normalization coefficient. Let us introduce a fractional parentage 
expansion of the HCH ^ onto the complete set of HHs YKa-yaiPc) foi' the A — 1 nucleons: 

yK-fip)= {AK-/\A-lK,j,J) Yk^^SPc) 

x{L,MLjm\LMLKScMs.^a\SMs)iT,MT^^T\TMT) (30) 

where {AK-flA — life 7c, is the fractional expansion coefficient. Multiplying Eq. ([30|l by YK^-ySPc)i 
^K^nimaTiO, ^A-i) ^^'^ three relevant Clebsch-Gordan coefficients and then integrating over angular variables and 
making summations over all projections of angular momenta and spins, we get 



UK-f - VA {AKj\A - li^eTc, 0- (31) 

Thus, the normalization coefficient Afxy is directly related to the fractional expansion coefficient into the core state 
on which the HCH is originally built. The knowledge of this relation is very useful for calculating the overlap integrals 
between the wave functions of nuclei A and A — 1. 

The overlap I^^l that enters Eq. ((29|) can be calculated using the fact that for an arbitrary function = 
0Qj^(l, ...,A — l)ipa2{A), where (f) is an antisymmetric fmiction, 

{AiJo,'\A^^) = (1, ...,A- l)^„,JA)\A^Po.) 

= S„^, -{A- l)((/.„; (1, ^ - 1)(^„^ (^)|0o, (1, - 2, A)(^„, {A - 1)), (32) 

where A is determined by ([5]) . The exchange term in this expression can be calculated using the fractional parentage 
expansion of the translation- invariant oscillator wave function oi A—1 nucleons. The details of these calculations and 
final analytical expressions are given in Appendix B. 



B. Matrix elements of central two-body NN interactions 



In this paper, only central two-body NN forces are considered, 

A 

V = Y. J2 v!f'\\r.-r,\)hti^,J). (33) 

i<Cj s,t=Q,l 

Here Pst{i,j) is the projector into states with spin s and isospin t of the pair of nucleons (i,j). To calculate the 
matrix elements of this potential in the HCH basis using Eq. (p6)) . one needs to know the matrix elements V^,^ in 
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the oscillator shell model cluster basis, given by Eq. ([28]) with O = V. The latter can be found using the fact that 
for an arbitrary function tpa = (p^i (1, ^ — ^)Va2 (^) that contains a product of antisymmetric function (j)^ 

-(A-1)(0„,(1,...,A-1)^„-(A)|F|0„,(1,...,A-2,A)^„,(A-1)). (34) 

Eq. p4|) contains direct and exchange terms that can be evaluated with the fractional parentage expansion technique. 
The derivation of formulae for V^,^ is given in Appendix C and the final expression for the potential V^^i^ is 

Kfa'/= E ^^^^^C^o^o.* (^"i'oll^^^llV'no^o), (35) 

where 

POO 

is the matrix element of the NN potential between the two-body oscillator wave functions i'noioi''')- The analytical 
expressions for coefficients C^r^^/g^j are given in the Appendix C. 

Finally, the hyperradial potentials VK-i.K''y' (p) are calculated by using Eq. ([26]) and applying the Laplace transform 
(P7)l to the shell model NN matrix elements V^^f . This gives 

V^,,K'y{p)-iMK'yUK,r' E E<n.iostipWr^'Ztv (37) 

where 
and 

"^V^nolQSt ^ 2 E E '^"^"oiost ^n'K'J' B^KJ- (39) 
K=0 k'=Q 



VII. APPLICATION TO ^HE 

The hyperspherical cluster model (HCM) proposed here, is tested for the simplest cluster nucleus, ^He = ^He + 
n. Only one, the lowest order, hyperspherical harmonics i^ooo(Pc) is retained in the HCH basis states to describe the 
core ^He so that the chosen hyperspherical cluster basis 

yKjip) = J^kI AiYooa{pc)VK.nlmariO,^A^l)), (40) 

where n = {K — l)/2, is orthogonal. The norm of this HCH and the coefficients C'^^^^i^^^ required for calculations of 
the hyperradial potentials are derived using the shell model technique from Appendix C that gives 

and 

^kk' -qh^9-4'i-3W X X, X QX fx ^0000,00^0 . ^OOOO.OOnoO \ 

'-n'^nolost - J (i + ^ ) O^^'dnofiOn'^^Qdig^o " Mi^fi ^Ono,0 o,0,0,k1,OOk'1,1 + ''«i,0 Jo,0,0,k'1.00k14 J 

^ (2.s + l)(2t + l) ^^^^^^ l|^|no?oiVA : 1)(00k'1 : l\j\n'aloNA : 1) (42) 
JVA ^ ^ 

To solve the hyperradial coupled equations, the computer code STURMXX has been used. 
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FIG. 1: Binding energy (a) and the r.m.s. radius (b) of ^He calculated in the HCM for different cut-off hypermoment Kmax- 
The HSFM calculations, performed in a restricted basis that contains only irreps X = {2i + 1, 0), where i = 0,1,..., (/5'max — l)/2 
(open squares), and in full basis (filles circles) are shown as well. The binding energy and the r.m.s. radius of the *He core are 
represented by dotted horizontal lines. 



The calculations are performed with the Volkov VI potential |16[ neglecting the Coulomb interaction. The standard 
value of the Majorana parameter m = 0.6 used with VI gives unbound ^He. For unbound states, convergence in total 
energy may not be achieved and the radius of an unbound nucleus should diverge. Therefore, keeping in mind future 
applications for heavier loosely-bound nuclei, the HCM calculations are also performed for m = 0.3. With this value 
of ni the NN interaction in odd partial waves becomes attractive leading to a bound ^He but keeping the separation 
energy of the valence neutron relatively small, at about 2 MeV. The resulting binding energies and r.m.s. radii for 
both m = 0.6 and m — 0.3 are shown in Fig. la and Fig. lb as a function of the cut-off hypermoment iiTmax — 2nmax + 1- 
The binding energy and the r.m.s. radius of the core '*He, calculated in the lowest approximation of the HSFM, are 
shown in these figures by horizontal lines (these values do not depend on to). 

One can see that with only twelve basis functions the binding energy and radius of a bound ^He corresponding to 
m= 0.3 has practically converged giving E = -30.66 MeV and (r^)^/^ = 1.96 fm. For m = 0.6, the binding energy 
still slowly decreases while the r.m.s. radius diverges as expected. For Xmax < 9, the HCM calculations are compared 
to the HSFM calculations performed within the full HH basis using the technique from Ref. [l^. For iCmax = 9, the 
difference between the HCM and the HSFM is within 1.5 MeV but the HSFM calculations involve 117 bases states 
as compared to only five states in the HCM. 

It has been shown in [l^ that, from the point of view of the kinematic rotation group Oa-i, clustering leads to a 
particular symmetry in the total wave function. According to flT\, if each of two clusters have a fixed internal Oa-i 
symmetry the expansion of the total wave function should contain only those HHs that correspond to the irreducible 
representations (irreps) (Ai+2i, A2, A3) of the orthogonal group Oa-i, where z = 0, 1, 2, 00. For ^He, the HSFM 
calculations of Ref. [iSj that include only irreps \ = {2i + 1, 0) give the binding energy that differs only by about 0.4 
MeV from the exact solution. Such calculations correspond to taking into account all possible cluster partitions in 
^He and they are also shown in Fig. la and Fig. lb. The difference between the binding energies calculated within the 
HCM and the HSFM with A = (2i -I- 1, 0) shows that other cluster partitions than ''^He -|- n and/or more complicated 
structure of the "^He core contribute about 1 MeV to the total binding energy of ^He. On the other hand, good 
agreement between the ^He radii calculated in HCM and HSFM indicates that the long-range behaviour is properly 
represented in the HCH basis. 

To investigate the long-range behaviour in more detail, the radial overlap integral Iij (r) between the wave functions 
of bound ^He and ^He has been calculated for ni = 0.3. The definition of this overlap, 

^■W = I](^mia|jTO,) / drYi*^{r) x{,i,{'f'"%Xi,X2,X3)\'i/]^':,{xi,X2,X3,r)), (43) 

includes the wave functions of '^He and ^He written in non- normalised Jacobi coordinates Xi, the last Jacobi coordinate, 
X3 = r, being the distance between the valence nucleon and the center-of-mass of ''He. Such a definition is consistent 
with the one used in different reaction theories and in microscopic cluster models. The derivation of the final expression 
for this overlap is given in Appendix D. Since the spin-orbital interaction in these calculations is not present, the total 
angular momentum j is ommited below. 
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FIG. 2: Overlap integral (a) for ^He, shown in linear (main graph) and logarithmic (inset) scales, and its ratio to the function 
—ihihf'^ {inr) (b) calculated in the HCM with VI and m — 0.3 for different cut-off hypermoment ifmax. The function ~-iKh[^^ (tKr) 
is shown in the inset of (a). 



The calculated overlap /;(r) is shovifn is Fig. 2a both in linear and logarithmic scales as a function of the cut-off 
hypermoment iiTmax- With twelve HCH functions the overlap Ii{r) practically converges for r < 6 fm. The r.m.s. 
radius of this overlap, shown in Fig. 3, has almost converged, giving 3.82 fm. At large r, the /;(r) should behave as 

Vbliir) ~iCiKh[^\iKr) r oo, (44) 

where k = •\/2/xe/?i, /x is the reduced mass of ^He -|- n, e in the neutron separation energy, equal to the binding 
energy difference in ^He and ^He, h'"^^ is the Hankel's function of the first kind and Ci is the asymptotic normalization 
coefficient (ANC). The function —iKh^^\iKr) is shown in the inset of Fig. 2a on a logarithmic scale and one can see that 

with increasing -fCmax the bevaviour of the overlap /; (r) approaches the trend given by (j44p . The ratio ili (r) / nh^^^ [inr) 
is shown in Fig. 2b for r < 12 fm as a function of -ftTmax- Although this ratio still deviates from constant value for 
r > 7 fm, the convergence achived at 5 < r < 7 fm with twelve HCH functions makes the unambigous determination 
of Ci possible. The C; value obtained at these r is 0.72 fm~^/^. The spectroscopic factor for the overlap Ii{r), defined 
as 

/•OO 

5 = 5 / drr^lfir), (45) 

is shown in Fig. 3. It is almost independent of i^max and equal to 1.2 which is slightly lower than the value of 1.25 
obtained in translation-invariant oscillator Ohcu shell model. 

The expansion of the ^He wave function onto the HCH basis (|40p should provide very similar results to the traditional 
microscopic cluster model 

^ = He ^ ^(^^Hc ^ g,^r)YUr)X^.^^r) (46) 

with the closed Os-shell oscillator wave function for ^Hc. In the present paper, the total binding energy of ^He has 
been calculated in the MCM as well for different oscillator radii b and the results are compared to those calculated in 
the HCM in Fig. 4 for both values of the Majorana parameter, m = 0.6 and m = 0.3. For m = 0.3, the MCM binding 
energy of ^He taken at its minumum, E = -29.99 MeV is about 0.7 MeV higher than the HCM result E = -30.66 
MeV. This is because the *He energy obtained in the lowest order approximation of the HSFM, E = —28.58 MeV is 
by the same amount lower than the expectation value E = —27.89 MeV for the ''He energy in the oscillator Os shell 
model basis. The neutron separation energy in both cases is practically the same, 2.08 MeV in HCM and 2.11 MeV in 
the MCM. The overlap integrals calculated in both models are also practically the same. The spectroscopic factor is 
1.2 is both cases and the ANC is 0.73 fm'^/^ foj. mCM as compared to 0.72 fm^^/^ Jq^. hCM. The calculated functions 
gi{r) are also very similar. The gi{r) has asymptotic behaviour given by with the same ANC, however it reaches 
the asymptotic form at significantly larger radii thus being less convinient to use for the ANC determination. For 
TO = 0.6, the results of the MCM and HCM cannot be compared directly because of the different boundary conditions 
imposed for unbound states in these models. In Fig. 4 the resonance energy is presented for the MCM, while the 
energy obtained with Ky^ax — 23 is plotted for the HCM calculations. The HCM result is lower. 
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FIG. 3: The r.m.s. radius and the spectroscopic factor of the overlap integral Ii{r) for ®He calculated in the HCM for different 
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FIG. 4: Binding energy of ^He calculated in the MCM for various oscillator radii b in comparison with HCM values. 



VIII. SUMMARY AND CONCLUSIONS 



This work has demonstrated that the HH basis can be reorganised in such a way as to represent the cluster 
structure of weakly-bound nuclei and the long-range behaviour of their valence nucleons. This can be achieved by 
antisymmetrizing the product of completely antisymmetric HH for the core and the relative hyperangular function 
of the last nuclcon. The hypcrsphcrical cluster model that retains only a few HHs for the core while performing 
summation over all quantum numbers related to the valence nucleon strongly resembles traditional microscopic cluster 
models. However, in the HCM the factorization of the wave function into the core and the relative wave functions 
is absent which reflects the influence of the last nucleon on the core wave function. Such influence is absent in 
microscopic cluster models where the wave function of the core docs not depend on the position of the last nuclcon. 
The factorization into core and relative wave functions can be achieved in HCM by projection onto eigenfunctions of 
the A — 1-body Hamiltonian. This makes the HCM to be equivalent to a multichannel microscopic cluster model with 
more complicated core structure and with all possible monopole core excitations built on HHs retained in the HCH 
expansion. On the other hand, the HCM is a narrower version of the HSFM because its hyperradial part is found by 
solving the same coupled set of differential equations that arises in the traditional HSFM. 

The feasibility of the HCM calculations has been studied using ^He as an example. With only one HH for the core 
^He retained in the HCH expansion, such a model should be very close to the microscopic cluster model with the 
core described by the lowest approximation of the HSFM because all the monopole hyperradial excitations lie too 
high in energy. The calculations, performed with a simple soft-core VI potential, have confirmed that it is possible 
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to get a converged solution for binding energy, r.m.s. radius, overlap integral and ANC with a small number of basis 
functions if the strength of the Majorana force is modified to bound '^He. The values obtained of these quantities are 
very close to those calculated in the traditional MCM with the oscilator radius that gives the minimum energy in ^He, 
as expected. For unbound ^He, the convergence of the HCH expansion can not be achieved and the r.m.s. diverges 
because of the wrong boundary conditions imposed on the HCM wave function. Scattering problem should be solved 
in this case rather than a bound state problem. 

The long-range behaviour of the last (bound) nucleon, governed by the Hankel function, is well reproduced when 
the number of HCHs increases. In the present short study, twelve HCH basis states are enough to derive the ANC 
although more functions are needed to reproduce the exponential decrease at r > 7 exactly. However, a further 
increase of basis states in the framework of HCM may not be the best option. The HCH basis can be used in the 
hyperspherical interpolation approach of Refs. and [TH where asymptotic conditions are incorporated into the 
system of radial equations, making the size of the basis smaller. Scattering problem can be considered in this basis 
as well. 

In the present test study, only the simplest core structure has been considered. Application to a case with a more 
complicated core, for example, the one with particle- hole excitations, can be done straightforwardly using formulae 
from Sec. VI and Appendix C. Currently, the only missing input quantities for such calculations are the FPCs. 
However, they can either derived using the recursive procedure of Ref. ; IQ] or by overlaping the HHs obtained in the 
shell model approach developed in Ref. [14]. 

Finally, the ideas presented in this paper can be extended to nuclei with two loosely-bound valence nucleons, and 
in particular to Borromean nuclei. This would include more complicated matrix elements and fractional parentage 
expansions but will make it possible to achieve a proper three-body description of such systems at large distances 
within a many-body formalism. 
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IX. APPENDIX 



A. Expansion of Jacobi polynomials into oscillator wave functions 



The Jacobi polynomial P"'^ {cos 29) can be represented as follows: 
On the other hand, 

K = ^ ' / V 

where 

= '"'^\/f55rTTT|y (!)' (^v^') (49) 

is the radial oscillator wave function with an arbitrary oscillator radius b and L^^t^^^ is the generalised Laguerre 
polynomial. Eq. can be proven by induction. Inserting and (I49p into (|T7)) gives relation between Jacobi 

polynomials and the oscillator wave functions. 



B. Normalization of hyperspherical cluster harmonics 

To calculate the exchange term in Eq. ((32|l the fractional parentage expansion of the oscillator wave function of 
the core can be used both in bra- and ket-vectors. This reads: 

x{L''M'lhmi\L,ML^){S''M'^^a,\S,MsJ{T''M^^T,\T,^^^ 

(50) 

where \A — 2N" K"^" M'lM'gMlf) is a translation-invariant shell model wave function for ^ — 2 nucleons with total 
number of oscillator quanta equal to N'\ well-defined hypermoment K" and other quantum numbers denoted by 
7". The fractional parentage coefficients (FPCs) {A — \NcKc"ic\A — 2N" K"'y" , Kih) for this expansion have been 
introduced in Ref. [20]. In Eq. ([50)1 Nc is the total number of oscillator quantum and Nc = Kc = K" + 2ki + li. 
These FPCs can be calculated using the ideas of Ref. [l^. Alternatively, they can be derived by explicit overlaping 
shell model wave functions with well-defined hypermoment, constructed in Ref. [U, for A — 1 and A — 2 nucleons. 
It could even be more practical to expand the wave function ^qk^^^^"^'^" (^i, €^-2) i^^to wave functions for A — 2 
without a well-defined value of hypermoment. Therefore, below, the index K" in the wave functions for A — 2 nucleons 
in all FPCs will be omitted. The index Nc will be omitted as well because Nc = Kc- 

The wave function of the two last nucleons in the bra-vector of the exchange term (j32p is expressed in coordinates 
^i^d while in the ket-vertor it is expressed in ^^_2 ^^'^ ^'a-i^ 

= (E + r.) - (51) 

Using the Talmi-Moshinsky transformation 



IV'Kiii(lA^2)V^K/(^A-i) : ioA/o) = E (^''''^i^i • Lo\A{A-2)\K,h^l : L„)\iJ.'A$'A-i)Mi'M'A^2) ■ ^oMq) (52) 
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where {k'1'k[1[ : Lo\A{A — 2)\kiIikI : L^} is the Talmi-Moshinsky coefRcient (the ordering and meaning of symbols in 
this coefficient is the same as in Ref. [2l[) and 2ki + li + 2k + I — 2k'i + l[+ 2k' + I' , we get for the overlap the 
following expression: 

N"f" 

X (^_^s+s"+s.+s:+T+T"+T^+T^ J2 {A- 1K,jM - 2N"j", K^h) {A - IKj'M - 2N"j'\ k[1[) 
xY,{-y'^^'''^"U{L''liLl;LcLo)U{L''l[Ll';L'^Lo){K'l'K[l[ : Lo\A{A ^ 2)\kiIikI : Lq) 

Lo 

(53) 

where U{jij2Jjz\ji232s.) = ji2j23W^(jij2 Jjs; ji2j23), 3 = V2j + 1 and W is the Racah coefficient. 



C. Matrix elements of central two-body NN interactions 

The matrix elements (y^Va' I^I-^V'q) between the antisymmetrized oscillator cluster shell model wave functions 
contain direct and exchange terms, 

{A'^^,\V\AiP^) = Vd„ -{A- 1)T4,. (54) 

The direct term 

Vd„- = (0a; (1, -.A- 1)^^ {A)\V\cj>^, (1, A - l)^^, {A)) (55) 
can be separated into two parts, 

A-l 

Vd^r = + = <5„„„, (0„; (1, ^ - 1)1 5^ V^Mo.^ (1, -.A- 1)) 

A-l 

+ (0„; (1, ^ - 1)^,^ {A)\ V,A\(t>c., (1, A - l)V'a. (A)), (56) 

i=l 

the first of which is the expectation value of the NN potential for A— 1 nucleons of the core in the standard translation- 
invariant shell model basis, 

v!j^}^5..Ja'&...A%!tiy'^'{(^^^ (57) 

and the second term Vj^J. being the folding potential between the last nucleon and the core. The exchange potential, 

= {K-Sl,-,A)^^.^{A-l)\V\cj)^,{l,...,A-2,A-l),^^,{A)) (58) 
can be represented by four terms, 

= + + + (59) 

according to the following separation of the two-body interaction potential: 

A A-2 A-2 A-2 

v^j = va-i.a + Y.^^^ + Y. ^'^-1 + ■ c^o) 

i<j 2—1 2—1 i<j 

The terms vjj^^ and Vex' resemble each other. They can be easily calculated by separating the wave function of the 
{A ~ \,A) pair in the cluster wave function This can be achieved by using the fractional parentage expansion 



50| of the core wave function ^'q^"^ " ^''(^i, ■■■t$,a-2) combined with the Talmi-Moshinsky transformation: 

~A 



^ — 2 

|-0Ki/i(lA-2)V'Ki(^yi-i) : LqMq) = 2_, (ki^iKi^i : ^o| ^ |»-o^oA^A : Lq) \il}nouXCi)'^NK{C2) ■ LqMq), (61) 

noloNA 
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in which 



^ ^< ^ ^ / 2(A-2) / 1 \^ rA-x+rA \ 



Then the wave function ' reads as follows 

^%\,..iLSTM,MsM^ A-i) = E 1^ - 2iV"7"Mi'M^'M^) - lif,7c| A - 2iVV, 

iV " 7 " K 1 / 1 M^' M 

^ j4 — 2 

no'oJVALo 

X E Xs„Ms„ToM.„ - 1, A) U{S"\s\;S,So)U{T"\t\-T,To) 
X E {L"M'ILoMlo\LMl){S"M'^SoMs,\SMs){T"M:^ToMt,\TMt). 

(63) 

The expansion of the wave function ^^^) in the bra vector looks exactly the same apart from an additional phase 
factor (— )'o+'So+To ^hich appears due to interchange of ^ — 1-th and ^-th nucleons. As the result, we get: 

Vi^-{A~1)V}^^ = {A-1){^'-^}{1,...,A-1,A)-<P'-^}{1,...,A,A-1)\VA^^^ 

rtQUolost 



where 



'^U{T"-T-; T,t) U{S"-S-; S',s) U{T"-T-; T',t) E U{L"hLl- L,L^) U{L"l[Ll'; L'.Lo) 

A — 2 A — 2 

Y,{kiIi,kI : Lo\^-\nolo,NA : Lq) {k[1[,k'1' : Lo\—-\n'olo,NA : Lq) 



X 

NA 



(65) 



K2) 
I' ex 



To calculate the exchange term Ve^ 

A~2 



K(?^ = (</><(l,-,A)(^a^(^- 1)1 E V^M|</>a,(l,...,^- 2,^- l)^„,(^)) 

i=l 

= (A - 2)((.^„, (1, A - 2, A)(p„, {A - 1)\Va-2,a\^c.^ (1. A- 2, A- l)<fiaM)) (66) 

it is convenient to separate the wave function of the pair [A — 2, A — 1) from the cluster function ■* using a 
two-nucleon frantional parentage expansion for the core, 

'""^"^ (€i: = E (^ - l^'-^-l^ - 3iV"/^"7"; 7VA, noZo^oTo(Lo) : Le^eT,) 

jV" l<-"7" JVAno io -C-o SoTo 

X E 1^ - ^N"K"L"M'lS"M'^r'M^)\^ljNA{V2),^noio{Vi) ■■ LoMlo){L"M'ILoMlJL,MlJ 

M'l M'^ M'^ Mlo Msq Mtq 

x{S"M'^SoMs,\S,MsJ{T"M;fToMT,\T,MTjxSoMs,ToM^,{A - 2, A), (67) 

where \A - "iN" K" L" M'lS" M'^T" M^) is the translation- invariant shell model wave function for A - 3 nu- 
cleons with N" quanta, hypermoment K" and other quantum numbers denoted by 7" and {A — lKcjc\A — 
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3N"K"j";NA,noloSoTQ{Lo) : LcScTc) is the two-nucleon FPC. As explained in the previous section, K" is omitted 
below. In Eq. ([S7)) . the coordinates r]^ and J72 are the following, 



A~3 



1 , . /2(A-3)/ 1 rA-2+rA^i \ 
r/i = ^(r^-.-r-^-i), = ^ ^-j—^Y.J^ ^ j. (68) 

Exchange of nucleons A—1 and A leads to new coordinates, ri[ and rjj, that are obtained from t]^ and J72 by replacing 
r^_i by rA- Applying the expansion (|67p and using the transformation 

^ {L''M'lLoML„\LcMLj{L,MLjm\LM)\i;NA{V2), i^nokXVi) ■ LqMlo)^^iU^a-i) = 

^ (L"M;:L;,'A/i'J4'A/i'J(4'AfIV"m"|iM) 

L''M'' L''M'' K"l"m"N'A'n'l' 

the wave function kILSTM^MsMt A — 2, A — 1, A) can be rewritten as 

^^kI^^ilstm^MsmJ^ ...,A-2,A-1,A)^ " ^ N'^'M'^M'^M^) 

N"'y" NAN' A'nolon'gl'gn"l" 
LoL'gL"SoS'Q S'jToTgT" {n} 

{A-lK,j,\A~3N"y';NA,noloSoTo{Lo) : L.ScT,) \^N'A'iV2)^n'xiv[) ■ L'oM'^„)ip^,,l''n^''{^'A-l) 

\S \ S, 

xxs^M', t;,m^ (A - 2, A)x.'r'iA - 1) S^S'^SoS'^XX'ToU { ^ I S'^ 

Sc So S" J I f , To T" 




(70) 



where {fi} = {M'^M'^Mf/^MLM'sM^Ms^MToML^miaT} and 

NAninihCC 

A_o ( L I" L'n ( L lo Xc ] ( L l„ X, ] 

x{n'^l'^N'A' :L'^\-—\nlnih:L'i)y^{X',t'K^^^ h L'^ ) { I A C ) { I AC W71) 

^" ^ A.A^ [ A', C L" J [ £' L" J [ io i" J 

Then in the wave function ^'^K'-y' n'v lstMlMsMt^^^ ■■■^ ~ 2,71, A — 1), in which the nucleons A—1 and A are 
interchanged, it is suffucient to extract the pair {A — 2, A). This gives 

'i'^!.^«'rLSTM.M.M. (1, ^ - 2, A, A - 1) = E 1^-3 WMi'M^'M^) 

Ar"7"JV'A'<i^'Li'S;,'T^'{A<'} 

x(A- lif^7cl^- 3iV"7";^'A',</[,'5^'T^'(L^,') : L'.S'X) l^^'A' (»72)^<;i' (rj'i) : (^:4-i) 
y-XS'^M'^Tl^M'^^^ {A ~ 2, yl)x.r(A - l)(L"A/i'L;;M£'JL^M£J(L^MiJ"m"|LMi) 

X I S',M's:}{S',M's^ \a\SMs){T" M'^T^; M'^^ \T'^M!r^ ) iT^M!r^K\T Mr) , (72) 

where = {M'lM'^.M:^M'l^M'^.^M^^^M'^^M'g^M:^^m''aT}. Using Eqs. dZQ]) and ^ in ^ gives the exchange 

(2") (3) 

term Vex ■ The exchange term Vex , 

A-2 

V}^^ = iK', (1, ^ - 2, A)¥'^ (^ - 1)1 E ^^^-1 l<^"i (1, ^ - 2, ^ - l)^„, (A)) 

1=1 

= (A - 2)(^„- (1, .4)¥>^ (A - l)|yA-2,A|</'ai (1, ...,A-2,A- 1)^„, (A))) (73) 
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can be calculated in a similar way to provide the final result 

Ki'^ + Ki'^ = - 2) E ^^i^Mst + W^nt<;i«t) {^n'il'^\VXU,A\\'^n',lO^ (7^) 

where 

W:^fn',i;st = E <^ - ^Kclc\A - 37V"7"; iVA, no^o5oro(Lo) : L,S,T,) 

N"j"N\N'\'noloLoL'aSoSl^ToT^ 

x{A- IK'.j'JA - 3N'Y; N'A', n'Xst{L'^) : L'XK) 

><7^™^i;4'^'^^ ^^^^^^ i ^ \{ \ i ' • (75) 

" \S, So S" ] {t. To T" ] 

Finally, the exchange term 

A-2 

Ki'^ = (<^a;(l,-,AK^(^- 1)1 ^y,,|0„,(l,...,A-2,^- l),<p„,(^)) (76) 

can be caclulated by separating nuclcon ^ — 1 in the bra and nucleon A in the ket vector using one-nucleon fractional 
parentage expansion and the Talmi-Moshinsky technique. This gives 

KL'^ = E K"l>.5s.s-^T"fM - 2A^"7"II E Vi,\\A - 2N"i') U{S"\\s; S,S',) U{T"\\t; TX) 

JV"7"JV"7" »<J 

X (_)S+s"+Se+s;+T+T"+T,+T^ E (^ - l^c7c| A - ^N'^ , mh) {A - IK'^-y'dA - 2N"'y", k[1[) 
X Y,i-y'^^'''^°U{L''liLl-LeLo)U{L''l[Ll';L'^Lo){K'l'K[l[ : Lo\A{A - 2)\kiIikI : Lq). 

Lo 

(77) 

This expression contains the shell model matrix elements of the potential energy in the A — 2 core, {A — 
2N"^"\\J2t<j^ Vij\\A - 2N"j"), that can be calculated in the usual way. 



D. Overlap integral 

The overlap integral /(r) = that enters reaction theories like the distorted wave Born approximation, 

coupled reaction channels, breakup and capture theories, is a function of distance between A'^ and the center-of-mass 
of ^ — 1. It is defined in terms of wave functions ^a-i and ^'a that depend on non- normalised Jacobi coordinates 
Xi = l/«X]}=i ''i " ''i+i- On the contrary, the wave functions in the hyperspherical formalisms can be rewritten in 
normalised Jacobi coordinates = aiXi, where = + 1). Since 

A-l 

dpA = d(,^...d£,A_^ = ]J a;^ dxi, (78) 

i=l 

the wave function ^'^({a;,:}) normalized in coordinates {xi} is related to the wave function ^A({^i}) normalized in 
coordinates {^j} in the following way: 

/A-i \ 1/2 

^A{{xi}) = ( n an ^A{{aiXi}) = A-^I^^A{{aiX,}). (79) 

Therefore, the overlap integral I{xa-i) used in reaction theories is related to the overlap integral I{^a_i) obtained 
in normalized Jacobi coordinates as follows: 

I{XA-l) = j dXi...dXA-2'^\{{Xi\)-^ A{{Xi}) 
= \-^) J d^,...d^A-l^\{^l,-,^A-2)^Ai$l,...,U-2,aA-lXA-l)=i^-^j I{^A-l)- (80) 
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Using the standard for reaction theories partial wave decomposition 

Ij[flTi^)= E (jm,,/A-iMA-i|JAAf^)(/miabm,)/,f'^--^(r)y,™(f)xi.i., (81) 

jnij Una 

where r = Ji{Mi) is the total momentum (its projection) of nucleus i, x^a^r ^^'^ spin-isospin function of 

the separated nucleon with spin (isospin) projection a (r), and expanding the ^He wave function in a hyperspherical 
cluster basis (PU]) . we get the final expression for the radial part of the overlap (*He(Xin|^He), 

+'i-=»-'>/ — ^ "--(Itf)- 

Here {K — 2n + l\K = Q ^l) is the FPC which can be calculated via the norm of the hyperspherical cluster harmonics. 
The overlap Ii{r) can be represented as a sum of direct and exchange terms. The direct term is an analog to the 
relative function gi{r) of the traditional microscopic cluster model Ij46[) and it can be obtained by removing the 
antisymmetrization operator from (|40p . The analytical expression for this term is obtained from (j82p by replacing 
the FPC {K ^2n + l\K = G,l) hy 5/(K = 2n + l\K ^ OJ). 



